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represent (M by M) positive definite matrixes, C an (M by N) matrix, D an (N 
by N) positive definite matrix, an inverse matrix of A is given by 

A-'=B-BC(I>C^BC)'c"B (33) 
Applying this Lemma to the calculation of the inverse matrix in Eq. (32), 
5 HL(n)HL"(n)+o^I=B-*+CD-'c", 

HL(n)HL"(n)<;D*'C", 

a^I=B-'. 

h(n)=C, 

D-'=I, 

S5 10 HL"(n)=C" 

jf By using them to calculate Eq. (33), the inverse matrix calculation in Eq- (32) 

J ^ is solved. Incidentally, Eq. (33) also contains an inverse matrix calculation 

}i (IHC^BBC)"^ but this inverse matrix is also scalar, and hence it can similarly 

! . be calculated. 

If 

|y 15 The thus calculated tap coefficient value G{n) is used to perfonn 

jS linear filtering of the difference value Rc(n) by the transversal filter 1 11 , and 
\^ the filtering result Z(n) and the esthnated channel impulse response value 

Bt(n) stored in the storage part 106B are used to calculate Eq. (25) in the soft 
decision value calculating part 120 to derive the soft decision value Ai[b(n)]. 
20 The equivalent amplitude n(n) in Eq. (25) is given by the following equation . ^ 1 0-^ 1 1 

^ysufetitutingEq:^^ EOji^^ 9(1 ^ V 

m = [HL(n)«[HL(n)HL(n)" WI]-^Hl(ii)] 
=HL(n)"G(n) (34) 
Incidentally, in the first iteration of equalization the tap coefficient calculating 
25 part 112 of the linear filter part 110 may also use the Matrix Inversion 
Lemma to calculate the tap coefficient G(n). 

As described above, this embodiment permits reduction of the 



number M neeck to be at least 2 or more. 

Assuming that the filter output can be approximated by a Gaussian 
distribution, Z(n) by Eq. (21) can be described as follows (see V. Poor and S. 
Verdu, 'Trobability of Error in MMSE Multiuser Detection," ffiEE Trans. 
Information Theory., vol. IT-43, No. 3, pp. 858-871, May 1997). 

Z(n)=M(n)b(n)+Tj(n) (22) 
where \i{n) is an equivalent amplitude of the ou^ut signal and Ti(n) represents 
a Gaussian distribution of an average 0 and the variance v^(n). Therefore, 
H(n) and v(n) can be expressed as follows: 
Si^^n)b(n)}^ 

,=lH,;;(n)«[H„,(ri)A(n^(n)«+g5^ 
where Qy indicates the element at the intersection of the J-th row and the J-th 
column of the matrix. 

v^=var {Z(n)}=M(n)-ji*(n) (24) 
From the above, the extrinsic infonnation derived by the linear equalizer can 
be deduced from the following equation. 

Xi[b(n)]=logP<^°>'^(") = -^^) = ^^^^HZ(n)} 

^p(Z(n)|b(n)=:-l) l-.n(n) ^^"^^ 

To obtain the optimum tap coefiScient G(n) with the above-described 
method, howevar, it is necessary that the inverse matrbc calculation expressed 
by the following equation be conducted in Eq. (20) for each point in time; 
hence, much time is required for calculation. 

dXji) = [H„,(n)A(n)H„,(n)« + o^I]-' (26) 
With the conventional method for iterative equalization by the linear 
equalizer, the inverse matrix calculation needs to be conducted for each point 
in time to update the tap coefficient-this gives rise to the problem of high 
computational complexity. 



